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Abstract 



We consider a two-component weakly interacting Bose-Einstein condensate 
in the presence of an external field which couples the two components. We 
express the Hamiltonian in terms of the energy eigenstates of the single-body 
part of the Hamiltonian. These eigenstates are the atomic dressed states of 
quantum optics. When the energy difference between the two dressed states 
is much larger than the mean-field interactions, two-body interactions in the 
dressed state basis that do not conserve the number of atoms in each of the 
two dressed states are highly suppressed. The two-body interactions then take 
on a simplified form in the dressed basis with effective coupling constants that 
depend on the intensity and frequency of the external field. This implies that 
the chemical potential as well as the quasiparticle spectrum may be controlled 
experimentally in a simple manner. We demonstrate this by showing that one 
may achieve significant variations in the speed of sound in the condensate, a 
quantity which has been measured experimentally. 



I. INTRODUCTION 

The recent experimental realization of trapped Bose-Einstein condensates (BEC) with 
internal degrees of freedom corresponding to different hyperfine states [0J || has sparked 
much theoretical and experimental study of the properties of multi-component condensates. 
Multi-component condensates exhibit a rich variety of new phenomena not present in con- 
densates with a scalar order parameter. For condensates with a vectorial order parameter, 
such as the F = 1 hyperfine multiplet of 23 Na, there have been predictions of spin waves, 
instability of vortices with more than one unit of circulation, coreless vortices, dynamic 
spin localization and the observation of spin domains 0. For two-component 
condensates, such as 87 Rb, there have been extensive theoretical and experimental investi- 
gations of the dynamics of the condensates in the presence of external fields which couple 
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the components. Examples include Ramsey fringes , non-linear Josephson oscillations || , 
instabilities in the quasiparticle spectrum [[J, the elimination of the mean-field shift in the 



energy [IC], control of the spatial dependence of the two components |jTT|, and collapse and 



revival of Rabi oscillations 12 



In this paper, a two-component homogenous weakly interacting BEC at zero temperature 
in the presence of an external field which couples the two states is investigated. The physical 
system is similar to that discussed in ||. In quantum optics, the dressed states, which are 
the energy eigenstates of a two level atom interacting with an external radiation field, form 
a convenient basis for many problems |13|]. More specifically, if one has a two-level system 
with states \a) and \b) governed by the Hamiltonian, 
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then the dressed states, |c) and \d), are simply the energy eigenstates of H with eigenvalues 
±|a; crf . The dressed states are a superposition of the states \a) and \b) with amplitudes 
determined by 5/Qr. We explore the use of a dressed state basis for describing an interacting 
two-component condensate. 

When the two-body interactions in the second-quantized Hamiltonian are rewritten in 
terms of a dressed state basis, one finds that there are terms that conserve the number of 
atoms in each of the two dressed states and terms that change the number of atoms in 
each of the dressed states. When the energy difference between the dressed states is much 
larger than the single particle kinetic energies and the mean-field energy, the Hamiltonian 
simplifies since the terms that change the number of atoms in each of the dressed states 
may be neglected. In this limit, the ground state of the condensate consists of atoms in one 
of the dressed states only, and the calculation of the excited states becomes trivial. In this 
case, there are two branches to the spectrum of elementary excitations. One branch has 
the standard Bogoliubov dispersion relation while the other branch corresponds to single 
particle excitations. In addition, the speed of sound is an explicit function of 5/VLr. This 
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result indicates that the spectrum of collective excitations of the condensate depend on the 
coherence between the internal states of the atoms in a manner which can be experimentally 
controlled. 

Before proceeding, it is helpful to review the hydrodynamics of a single component BEC 
at zero temperature. At zero temperature, the thermal component of the Bose gas is absent 
and the quantum depletion is negligible if the bosons are weakly interacting. In this case, 



the wave function (or order parameter) for the condensate, 0(r, i) = yn(r,t)e lS ( r,t \ obeys 
the time-dependent Gross-Pitaevskii equation [[14]] which may be written as two coupled 
equations for the density, n(r, t), and velocity, v(r,t) = — V5'(r ) t) 

d 

— n + V • (vn) = (1) 



and 



d / , N „ Ti 2 _ ~- mlvl 2 



m _ v + V ^(r) + I7n- s ^V»^ + -Li-j=0. (2) 

Here U = 4n ^ a , a, is the s-wave scattering length which characterizes the two-body interac- 
tions between atoms, m is the atomic mass, and V(r) is an external potential. In general, 
the quantum pressure, 2 J^ V 2 a/^, may be neglected in comparison to the mean-field inter- 
action. For a static condensate ground state, this corresponds to the Thomas-Fermi limit 
for the density, 



n, 



,(r) = U-\n-V(r)), (3) 



where \i is the chemical potential. By considering small fluctuations, 5n(r,t), about the 
static ground state density n (r), one may derive a linear wave equation for 8n(r, t) from 
Eqs. {§■§, 

d 2 5n(r, t) = V • [u 2 (r)V5n(r, t)] , (4) 

where u 2 (r) = ^n Q (r) is the local speed of sound. The sound speed may also be obtained 
from the relation, w 2 (r) = where P is the pressure of the ground state of the con- 

densate. The derivation of Eq. (f|) neglects the quantum pressure term in Eq. (||]). When 



V(r) = 0, the condensate ground state is spatially homogenous so that n G (r) = n . In 
this case, the qauntum pressure may be neglected when the variations in the density are 



small over distances on the order of the healing length, £ = 1/ ^y8nan , or equivalently, the 
elementary excitations have momenta satisfying hk <C mu where k is the wave vector for 
the excitation and u = u(r). 

For a homogenous BEC (i.e. V(r) = 0), the energies of the elementary excitations form a 
continuous spectrum. The dispersion relation for the long wavelength collective excitations 
of the condensate have the form u = uk <C These excitations correspond to phonons. 
The speed of sound in a homogenous BEC (sometimes called Bogoliubov or zero sound) was 



first derived by Bogoliubov |16| using a microscopic theory of weakly interacting bosons, 



u 



IAttK an , 

2 — • ( 5 ) 

m z 



The Bogoliubov theory is used in this paper to derive the excitation spectra for a two 
component condensate. 

Bogoliubov sound in a condensate of sodium atoms has been experimentally studied 
by M. R. Andrews and coworkers in a highly elongated cigar shaped trap jfT/j] . Note that 
V(r) 7^ in this case but along the long axis of the trap, the condensate can be treated 
as being locally homogenous. In this experiment, the repulsive dipole force of a focused 
blue-detuned laser beam was used to create a localized density perturbation at the center 
of the trap. The double peaked perturbation was subsequently imaged as it propagated 
along the long axis of the trap. When the density of atoms in the trap was varied, the 
measured sound speeds showed good agreement with Eq. fl5|) where n refers to the density 
at the center of the trap divided by 2. The extra factor of a 1/2 has a simple interpretation. 
Since the sound propogation is confined to the long axis of the trap, the condensate can be 
treated as an effective one-dimensional system. In this case the effective density is given by 
averaging the condensate density over the directions transverse to the long axis which in the 
Thomas-Fermi limit is just the central density divided by 2 WE . 



This experiment has been analyzed theoretically by several authors [ 18j [19| in the 



Thomas- Fermi limit. In a trapped condensate, sound waves can propagate provided that 
the excitations satisfy both hk <C mu(r) and fcl > 1 where L is the size of the condensate 
in the direction of propagation. The latter condition corresponds to wavelengths much less 
than the condensate size so that the variations in n (r) over distances on the order of the 
wavelength are negligible. Consequently, the condensate can be treated as locally homoge- 
nous. As such, the results derived for a homogenous BEC using Bogoliubov theory may be 
used for trapped condensates with the substitution n Q — > n (r) where the bar denotes the 
possible averaging over transverse dimensions for anisotropic condensates. 

In the following section a second-quantized Hamiltonian for a two-component BEC is 
given in the "bare" atomic basis and the Hamiltonian in the dressed basis is derived. In 
section III, the condensate ground state and elementary excitations in the dressed basis are 
determined. Finally, in section IV, the results are discussed including justification of the 
physical approximations used. 



II. PHYSICAL MODEL 

We consider a collection of N bosonic atoms that have internal states \a) and \b) with 
energies hu /2 and —Hu /2, respectively. There is a spatially uniform radiation field with 
frequency u e which couples the two internal states with a Rabi frequency Qr ||20|| . The 
atom field detuning is denoted by 5 = uj — u e . The atoms in states \a) and \b) may also 
be subject to an external trapping potential V a (r) and V&(r), respectively. Furthermore, the 
atoms interact via elastic two-body collisions through the interaction potentials V^(r — r') 
for i,j = {a, b}. The many-body Hamiltonian operator describing the system is given by, 



H — H at om + HcoU 

#,„,„„ = I d 3 rUl(r) 



(6a) 



2^ V + ^< r)+ T 



^(r) 



+hQ R [*l(r)* 6 (r) +*J(r)* a (r)]} 



(6b) 



H coll = ^J d 3 rd 3 r' {*t ( r )^rt ( r ')y ao ( r - r')^ a (r / )^„(r) + ^(r)^^)^^ - r / )^6(r')^(r) 
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+2& a (r)¥ b (r')V ab (r _ r ')# a ( r )^( r ')} . ( 6c ) 

Here, H atom is the single particle Hamiltonian and H co ii represents two-body collisions. 

The operators ^(r) and ^|(r) are bosonic annihilation and creation operators for 
an atom in state i = {a, b} at position r which satisfy the commutation relations 
[Vi(r), ^t(r')] = SijS (r — r') and [V$(r), ^-(r')] = 0. The operators, ^i(r), have been writ- 
ten in a field interaction representation which is rotating at the frequency of the external 
field, u e , 

y a (r) = ¥ a N \r)e 1 ^ 2 ; (7a) 
£ 6 ( r ) =*W( r )e- <w -*/ 2 ; (7b) 

where ^ (r) are the field operators in the normal representation. This explains the ap- 
pearance of the detuning in Eq. (|6bD . 

The two-body interaction between atoms, V^(r — r'), depends on the internal states of 
the atoms. For a dilute gas such that fi | a^- 1 3 <C 1 where n is the average density and is 
the s-wave scattering length between atoms in states i and j, the interaction may be written 
as Vij(r — r') = UijS (r — r') where C/y = 47r ^ a ' J . It is assumed that > corresponding 
to repulsive interactions. 

For the remainder of this section and sections III and IV, we consider only the case 
when Vi(r) = 0. One can expand the field operators in a basis of single particle momentum 
eigenstates, 

*.(*) = 4r£««n.e**/* (8a) 
v v p 

^( r ) = 7^E«^ ip - r/? \ (8b) 



p 



where V is the quantization volume and [a ip , a] /] = <%5 PlP ' and [cti P , (Xjp<) = 0. In this 
basis, the Hamiltonian becomes 

H a tom = £ | + y) a aP Qa P + ~~ Y ^ "Ip a bp + («fep"ap + u\pU b p) j , 

H C0 U — —rz 2_^i \^aa a api a ap2 a «P3 a ap4 ^bb a bp\ a bp2 a bp3 a bp4 + ^U ab a api Oi bp2 Oi a p 3 a fep4 

Pi+P2=P3+P4 



At this point it is advantageous to introduce the grand canonical Hamiltonian, K = 
H — fiN where \i is the chemical potential and N = J2 P { a \ P a a P + a lp a bp) is t ne number 
operator. The motivation for using K is that the Bogoliubov prescription which will be 
used later, results in a Hamiltonian that does not conserve the number of particles. The 
chemical potential serves as a Lagrange multiplier which allows one to impose the constraint 
(n^j = N ||21|| . On the other hand, since fi is the energy of an atom in the condensate ground 
state, K corresponds to a representation in which all single particle energies are measured 
relative to the condensate. 

Dressed operators c p and d p , which correspond to the atomic dressed states |c) = 
cos 9 | a) + sin (9 |&) and \d) = — sin^la) + cos9\b), may be introduced which diagonalize 
H a tom — I^N. They are related to the operators «b p and a ap by 



y dp j 



( 



V 



cos 9 sin 9 



sin v cos i 



\ 



/ 



a 



ap 



(10) 



^ &bp j 



1/2 



or tan 29 = 20a. 



where the dressed state angle 9 is given by cos 9 = -4= ( 1 H — , = 2 
Note that 9 is a measure of the relative mixing of the atomic states \a) and \b) in the dressed 
states so that, for example, 9 = tt/4 (Qr ^> \6\) corresponds to dressed states which are 
equal superpositions of the states \a) and |6) . In contrast, for 9 = (5 > and Qr ^1^1) 
one has \c) = \a) and for 9 = n/2 ( 5 < and Qr \8\) one has \c) = \b) so that in these 
limits, the dressed states may be identified with the atomic states. It is easy to show that 
[c p , c p ,] = [d p , d v ] = 5 p>p > and [c p , c p /] = [d p , d p >] = [d p , c pl ) = [d p , c p >) = 0. 
In the dressed basis H a tom — H-N has the form 



hu cd \ + / p 2 



cd 



fi d^dp \ ; 



(11) 



where hu> c d = hy5 2 + 4f2|. is the energy difference between the two dressed states. However, 
in the dressed basis H co ii has a significantly more complicated form, 

Hcou = Yl { u i c lA 2 c P3 c P4 + u 24>i4>2 d Ps d Pi + u 34>i d U c Ps d P4 

Pl+P2=P3+P4 



+^4 (4>i4>2 d P3 d P4 + 4i4 2 c p3 c p 4 ) + u 5 (c^c^c P3 rf P4 + C? P1 C L C P3 C P4) 

+C/ 6 (c pl rf p2 ^P3 rf P4 + 4>l4>2 rf P3 C P 4 ) } • (12) 

The collisional couplings C/j are given by, 

\J X = U aa cos 4 6 + C/ 66 sin 4 6 + -C/ a6 sin 2 26; (13a) 

U 2 = U aa sin 4 + C/ 66 cos 4 9 + ^C/ a6 sin 2 26; (13b) 

C/ 3 = (f/aa + U bb ) sin 2 2# + 2U ab cos 2 2£; (13c) 

C/ 4 = \{U aa + Efo - 2U ab ) sin 2 2#; (13d) 

C/ 5 = sin 29{U bb sin 2 - U aa cos 2 + U ab cos 20); (13e) 

C/e = sin2#(f/ 66 cos 2 # - U aa sin 2 6 - C/ a6 cos20). (13f) 

Even though the atom-radiation field interaction term has been eliminated, it would 
seem that nothing is gained by using the dressed basis because of the increased complexity 
of H C oii- However, one may make a transformation to an interaction representation with 
respect to the internal dressed state energies by defining new slowly varying operators c p 
and dp which are related to the normal dressed operators by, 

c P = c p e~^ 2 ; (14a) 
d p = d v e +l ^ dt/2 . (14b) 

The time evolution of c p and d p is governed by the Hamiltonian, K s , given by 

ks = £ { ~ i) 4>c P + - v) tfpdp} + 

^77 E {^iCpiCp.CpgCp, + C/ 2 J P1 J P2 J P3 J P4 + C/3C Pl rf P2 c P3 J P4 

Pl+P2=P3+P4 

+C/ 4 (e 2 — 4 c pi c P2 J P3 J P4 + e- 2 — *J P1 J P2 c P3 c P4 ) + £/ 5 (e^*c pl c p2 c P3 J P4 + e"*--* J P1 c p2 c P3 c P4 ) 
+C/ 6 (e^*c Pl J P2 J P3 J P4 + e-^* J Pl J P2 J P3 c P4 ) } . (15) 



The terms proportional to C/4, C/5, and oscillate with a frequency that is a multiple of oj C( i- 
The time scales over which the operators c p and d p evolve are governed by the single particle 
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2 

kinetic energy, Jp, and the mean-field interaction energies, n a Ui, which are on the order of 
noUaa, n Ubb, and n U ab where n Q is the density of atoms in the condensate. Consequently, 
when the dressed state energy splitting is sufficiently large so that 

P 2 

huj cd > — , n U aa , n U bb , n U ab , (16) 
2m 

the oscillatory terms in K s will undergo many oscillations before c p and d p will have changed 
significantly One may therefore consider the Heisenberg equation of motion for c p which is 
averaged over a time interval T = 2n/u c d, 

1 t-t+T dr 1 / >t + T 1 

fl -W dt =rl T^ dt - < 17 > 

This allows one to define a coarse-grained derivative, 

r^j, Mf + r)- epW (18) 

Jt at T 



5c p _ 1 ft+T dc r 
Jf~f 



When condition ( |i6|) is satisfied, the operators in K s may be treated as constant over the 
period T so that 

where K sr is the "resonant" Hamiltonian, 

k " = ? { \L ~ ") ^ + {h->) ^} + 

4? E {^iCpi^LVp* + U 2 d pi d p2 d P3 d Pi + t/3C Pl 4«2 5 P3^P4}- (20) 

Pl+P2=P3+P4 

Since c p changes very little in a time T one may identify the coarse grained derivative with 
the actual derivative, 

^£p ^ ddp 

The same results hold for d p . As a result, one may consider the time evolution of c p and d p 
to be governed by K sr . 
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Neglecting the oscillatory terms in K s is analogous to the rotating wave approximation in 
quantum optics |22[ p3j and has a simple interpretation in terms of energy conservation. As 
an example, consider the term cj^cj^dpg J P4 in K s . This term corresponds to the destruction 
of two atoms in the dressed state \d) and the creation of two atoms in state |c) which requires 
an amount of energy equal to 2hu c d. However when condition (|T6|) is satisfied, the kinetic 
energy and mean-field energy of the atoms is insufficient to overcome the energy difference of 
2fiui cd . Consequently, the process described by c ^ c l 2 d P3 d P4 can only occur over time scales 
consistent with the uncertainty relation AtAE « h, which in this case is At « l/2a> cc z. 
However, by coarse graining the time evolution of the operators c p and d p , one ignores 
processes that occur on time scales less than T > At. 

The slowly varying operators were useful for showing that several of the terms H co ii could 
be neglected. However, for the remainder of the paper, all calculations will be carried out 
using the operators c p and d p in the normal Heisenberg representation with the "resonant" 
Hamiltonian Kr given by, 

*- ?{£ + ^-'')M£- ! r-''M + 

7^7 £ {^ic Pl c P2 c P3 c P4 + f/ 2 d Pl 4 2 d P3 d P4 + f/ 3C Pl rf P2 c P3 rfp 4 }- (22) 

Pl+P2=P3+P4 

Notice that Eq.(|22|) resembles Eq. ( |9a}]9b"|) except that the external field coupling term is 
absent. Equation ( p^) is the central result of this paper. The coupling constants Ui, U 2 , and 
U 3 in Eq. ( p2|) are explicit functions of the dressed state angle. Consequently, the collective 
properties of the condensate which depend on the two-body interactions will be a function 
of 9. In the following section, the condensate ground state and elementary excitations above 
the ground states are calculated using Eq. (E3). 



III. CONDENSATE AND EXCITED STATES 



The c-number equations for the wave function of the condensate ground state may be 
derived by dividing the field operators ^ c ( r ) and ^d(r), in the dressed basis, into a part with 
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a non-zero expectation value and a fluctuating part, 5^(r), that has a vanishing expectation 
value with respect to the ground state, 

fc c (r) = -2° +<$* c (r); (23a) 

V 



* d (r) = -^ + ** d (r). (23b) 



The wave functions for the condensate ground state are <fi c = — ^ an d <fid = 

A pair of coupled equations for C , 0^, and /i may be derived from the expectation value 
of the Heisenberg equations of motion and using the fact that <fi c = <fi d = with respect to 



•d 



h l 

{J,<t>c = {^cd + Wd 2 + ^UMd\ 2 )(f>ci (24a) 



V4d = (-~^cd + U 2 \M 2 + lu 3 \<p c \ 2 )<i) d ; (24b) 



2 IT 1 2 

-u;^ + f/ 2 |0d| 2 + 

with the constraint n = ^ = |0 C | 2 + |0d| 2 . Equations ( |24a| - |24bD may also be derived by 



requiring that the energy of the ground state be an extremum with respect to variations in 
C and <j) d . 

Due to the non-zero mean-field interactions, there are three possible solutions to Eqs. 
(H|- gig ): 



>d 



2 h 

= n a , <fi d = 0, and /i c = -uo cd + lqn ; 

h 

2 = n , <fi c = 0, and fi d = ~^cd + U 2 n ; 



/...x i , ,2 {U 2 - ~Uz)n - Hud 2 (Ui - W 3 )n + hucd 

H M = c^K^ • W = + • 

1 (2EW 2 - 5^) n, + (Oi - t/i)ftw rf 
a ' ld " = 2 C/ 1+ C/ 2 -C/ 3 ' 

The derivation of solution (iii) requires <fi d ^ and <fi c ^ 0, but it is easy to show that in the 

limit that \(j) d \ 2 -> or |0 C | 2 -> one recovers solutions (i) or (ii), respectively. When ([II]) is 

satisfied, (iii) is an unphysical solution since either |0 C | 2 < or \<fi> d \ 2 < and consequently, 

case (iii) may be ignored. The solutions (i) and (ii) correspond to atoms of the condensate 

being in one of the dressed states, \c) or \d). Case (ii) corresponds to the thermodynamic 
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ground state since it has the lowest ground state energy. However, in what follows the 
excitation spectrum for both cases is calculated. 

The Bogoliubov approximation |16[] may applied to Kr for the two condensate solutions 
(i) and (ii) to obtain a Hamiltonian which is quadratic in the operators for the excited states 
(those states with p ^ 0). For case (i), one obtains the linearized Hamiltonian, 



K Rc = K oc + | + ^i^) 4> c p + \ U ^ n o ( c p c -p + 



CpC_p 



(25) 

by taking cj, ~ cq = y/No and d^ ~ do = and neglecting terms which are cubic and quartic 
in operators with p ^ 0. Here K oc = E oc — fi c N Q and E oc = ^f^No + ^UiNl is the total 
energy of the ground state. 

One can see directly from Eq. (^5|) that the excitations which are in the atomic state 
that is orthogonal to the ground state, i.e. \d), are single particle excitations with an energy 
spectrum of 

P 2 1 

huj d (p) = 2^~ huJ cd ~ {Ui - -U 3 )n . (26) 

This spectrum corresponds to the free atoms in state |<i) with a mean field shift of \U 3 n 
due to the interactions with the condensate. The term — {TiuJd + Uin ) in Ud(p) results from 
the fact that the single particle energies are measured relative to fi c . The reason c^d(p) < 
is due to the fact that ground state (ii) is the global minimum in the condensate energy. 

For excitations in state |c) , one can carry out a canonical transformation by defining 
quasiparticle operators, C p , such that 

J2 { (tt- + Uwo) c P Cp + \u x n (cpclp + c p c_p) 1 = ^ c (p)C P Cp + E vaCiC ; (27) 

p^O I \ Zm ) 1 J p^O 

where E vac ^ c is the vacuum energy for the quasiparticle vacuum. The C p obey the bosonic 
commutation relations [C p , Cy] = 5 PtP > and [C p , C p i\ = and can be expressed in terms c p 
and cp as 
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Cp = cosh (ppCp + sinh (p p cL p ; (28) 



The solutions for the quasiparticle energies and <p p are easily found to be pi 



P 2 ( P 2 



^{^ + 2U ^j- (29a) 

tanh2y? p = 2 UlUo — . (29b) 

2 

The long-wavelength quasiparticle excitations for which Jp <C 2U\n correspond to phonons 
with the dispersion relation hu c (p) = u c p and a speed of sound u c given by 



u r . 



'UWo 



rn 



(V aa cos 4 + C/ 66 sin 4 6 + hj ab sin 2 2fl) ; (30) 



The speed of sound given by Eq.(|30D is the same as that calculated from the pressure, 
P, of the condensate using P = along with u 2 c = ^J^-- It is easy to show that 

the excitations in state |c) give rise to density perturbations, 5h(r,t). The number density 
operator is 

n(r) = *t(r)* c (r) + ^t (r) ^ d(r) = _L £ ( C t /Cp + d \ pldp ) e^'^ h (31) 

p,p' 

which may be linearized around the condensate ground state to give n(r) « n + 5n(r, t). 
The density perturbations can be expressed in terms of the quasiparticles as 

5n(r,t) = (^) 1/2 ( cosh ^P - sinh ^ P ) (Ct(t)e-^ n + C p (t)e^' h ) , (32) 

V V J P7^0 

which has the form for the phonon states, p -C mu c , of 

1/2 

p«mu c \2Vmu c , 

Equation (|33| ) is the operator for the density perturbations which obeys Eq.(^) for a ho- 
mogenous fluid |21|. Consequently, the excitations in state \d) do not contribute (at least to 
order \fW ) to sound propagation. 

The calculation for case (ii) proceeds in an identical manner and, as such, we quote only 
the main results. The linearized Hamiltonian is, 
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1 /2 

Sh(r,t)= J2 ( y° P ) (C| ) (0)e- i ( p - r ^^/' i +C p (0)e ,(p ' r ^ c!,1)/R ). (33) 



K Rd = K od + | + U 2 n}j d p d p + ^U 2 n (d p dL p + d p d_ p ) j 

+ E (j^ + ^ - (^2 - ^sHI c i c p! ( 34 ) 



with if od = E od — yuA^ and -E od = —^f^-No + ^U 2 N% is the ground state energy. 
The excitations in state |c) are single particle excitations with an energy spectrum of 
hu c (p) = 2^ t + huj c d — {U2 — \U'i)n . The quasiparticle operators D p = cosh </? p <i p +sinh (fpdLp 
diagonalize the terms in KrcL involving dp and dp, 

E{g + d p d p + hj 2 n (rft dL p + V- p ) } = E ^(p)^p^p + (35) 

with quasiparticle energies 

hu d (p) 



P 2 1 P 2 



The long wavelength phonon excitations have a sound velocity u d given by 



«d = \/ = -\ — ^ sin 4 9 + U bb cos 4 + -U ab sin 2 20 . (37) 
\ m \ m \ 2 1 



Again, the sound speed given by Eq. ( |37| ) is identical to that calculated using P 
and u 2 d = t^J^-- The linearized density perturbations are now given by, 



dE, 



a,/ 



dV 



Sn(r,t) = ('^) V2 J2 (cosh^ p - sinh^ p ) (l)t (f) e -ip*/» + D p {t)e^ n 



1/2 

p<^mu d x - zv lnu d ' 

Equations ( |33"D and (|38[) for 5n(r, i) and Eqs. (|30D and (|37|) for the speed of sound are 



]T C n ° P ) (D p (0)e- l{p - r - Udpt)/h + D p (U)e i{ P- r - u ^/ n ) ; (38) 



the main results of this section. The expressions for u c and u d given by Eqs. ( p0[) and 
(ftT|) indicate that the speed at which a density perturbation propagates in the condensate, 
depends on the particular dressed state that the condensate atoms are in and the value of 
the dressed state angle. In the limit that U aa = U bb = U ab = 4w ^ a , one obtains u c = u d = u 
where u is given by Eq. (|5|). Consequently, the inequality of the scattering lengths is 
crucial for this effect to be observed. Figures l(a-c) show plots of u c and u d for various ratios 
of U bb /U aa and U ab /U aa . 
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The atoms in the condensate can be prepared in either states (i) or (ii) by adiabatically 
turning on the external field. Suppose that at t = — oo, = and 5 ^ so that the 
atoms are in either \a) or \b) . If the atoms initially in state \b) = \d) (5 > 0) , then the 
atoms will remain in \d) provided 9 varies sufficiently slowly For free atoms, the atoms 
will adiabatically remain in the initial dressed state provided l^cu^ 1 <C 1. However, for 
interacting atoms in a condensate, Ti\9\ should be mush less than all of the energies in the 
Heisenberg equations of motion which govern the time evolution of c p and d p . Consequently, 

2 

one must satisfy the more stringent requirement fnOd , n U aa , n Q Ubb, n U a b 3> ^|#|- 



IV. DISCUSSION 



Up to this point, the oscillatory terms in Eq. (p~5|) have been neglected under the as- 
sumption that their effect is small compared to the terms in K sr . However, it is important 
to obtain an estimate of the leading order correction due to these terms. The lowest order 
effect of the oscillatory terms is an energy shift in the chemical potential and quasiparticle 
spectrum which is analogous to the Bloch-Siegert frequency shift in the atomic resonance 
of a two- level atom due to the counter- rotating terms in the atom-field coupling [p^l . The 
shift in the chemical potential (i.e. /x — > fi + S/i where 8fi is the shift) for cases (i) and (ii) 
may be calculated by including the two-body interactions proportional to U{ for i = 4, 5, 6 



in Eqs. fl24aH24b|) . Doing so, one finds that for case (i) the shift is 8/i c = 3 ^ 5 J°- > and for 



-3(U 6 no) 2 



case (ii) the shift is SfXd - lh 

To obtain the shift in the quasiparticle spectrum, one may write Eq. fll5|) as K. 



s 



K sr + Hcs t nr(t) where K sr is given by Eq. fl2~0|) and H cs ,nr{t) consists of the terms which 
oscillate at uo c d- To calculate the energy shift, one first linearizes H CS)nr (t) around one 
of the ground states given by cases (i) or (ii) so that H CSynr (t) is quadratic in c p and 
dp for p ^ 0. Since e tuJcdt varies rapidly compared to c p and d p , the coupling be- 
tween the operators c p and d p in H CSjnr (t) may be adiabatically eliminated by integrat- 
ing the Heisenberg equations of motion for the quasiparticles and substituting the solu- 



16 



tion back into H cs ^ nr (t). One then finds for case (i), that Eq. ( plf) now has the form 



K Rc = K oc + E vac + hJ2 P ^o (w c (p) + Su c (p)) CjCp + (uj d (p) + 5uj d (p)) tf dp 



where the low- 



est order energy shifts of the excited states, including the shift in the chemical potential, 
are given by Mu c (p) = [ 7 |~,^ ) and Mw d (p) = (^n f-3(u s n f/2 ^ Note that the 



Ahuicd I huj c (p) j d,\fJ 2Tiu) 

p 2 

factor 7 ' 2 ^ j ^p 1 "" in Mu; c (p) comes from expressing the c p operators in the linearized H CS)Tir (t) 
in terms of C p using Eq. (|28|). For the energy shift to be negligible for the phonon states 
(i.e. lo c (p) 3> 5c<j c (p)), one must satisfy the condition ^ 3> ■ In a similar manner one 

may calculate the energy shift in the excited states corresponding to ground state (ii). In 
this case one finds that the energy shifts in Eq. fl34]) are hduj c (p) = — ^ 4 "°^ ^L^ 6 """* ~ anc ^ 
Ti5ud(p) = — ^12°] ^ 72 ^L d (p)"° ^ ■ Again, one finds that for case (ii), the energy shift for the 
phonon excitations is negligible provided ^ 3> . One can see that the energy shifts 

are smaller than the unperturbed energies by factors which goes like "C 1 and that even 
for the low energy phonon states, the energy shift is negligible. 

The validity of the results in the beginning of this section as well as in sections II and 
III all rely on condition ([H]) being satisfied. Since we are primarily interested in the low 

2 

energy phonon states for which <C n U aa ,n Ubb,n U a b, condition (|T6| ) reduces to 

hu cd > n U aa , n U bb , n U ab . (39) 

To satisfy this condition as the external field turns on, it is necessary to have h\S\ ^> 
n U aa , n Ubb, n U ab . At later times, however, once Qr is established, the detuning can be 
varied adiabatically to adjust the dressed state angle. For these times, to satisfy ( |3~9"1) it is 
sufficient that 

h2Q R > n U aa , n U bb , n U ab . (40) 
To estimate the mean-field interactions, one can take ~ a and use the values of a = 



2.75nm for 23 Na and a = 5.77nm for S7 Rb and a condensate density of n Q ~ 10 14 cm 3 [[15 



The mean-field energy, inh2an ° ; divided by h for 23 Na is then 9.5 x 10 3 s 1 and for 87 Rb it is 
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5.3 x 10 3 s _1 . Notice that this is comparable to the two-photon Rabi frequencies currently 
used in experiments with 87 Rb which are typically ~ 2n x 600s _1 |[20| . 

In order to obtain some feeling of the necessary field strengths needed to satisfy ([TED , 
one can, for the sake of definiteness, consider single photon rf transitions between Zeeman 
states in a spinor condensate so that hQ^ ~ ^bB where /ig is the Bohr magneton and 
B the magnitude of the rf magnetic field. Thus a magnetic field of B ^> 10 _3 (r would be 
required to satisfy (f40D . One could also hold 5 constant and vary the rf field strength B. The 
splitting between Zeeman states in the presence of a static uniform magnetic field B D will 
be on the order of uj ~ fisBo/h « 9 x lO 10 ^^ 1 for B D = 10 4 G. One could therefore achieve 
detunings such that 5 <C uo so that the rotating wave approximation remains valid but 
still have h \8\ 3> n U aa , n a Ubb, n U a b. Consequently condition (|39|) should be experimentally 
achievable. 

As mentioned before, in order for the speed of sound to show a dependence on the 
internal states of the atoms, the scattering lengths a,ij must be different. For magneti- 
cally trapped 87 Rb with the two 5SV/2 hyperfine states \a >= \F — l,m = — 1 > and 
\b >= \F = 2,m = 1 >, the scattering lengths are nearly equal and are in the ratios 
a-aa '■ a ab '■ a bb 1-03 : 1 : 0.97. This would preclude the use of 87 Rb. However, the recent 
demonstrations of the manipulation of the scattering length in condensates of 85 Rb and 23 Na 
using a Feshbach resonance p4| opens up the possibility of using similar techniques in mul- 



ticomponent condensates. An alternative method for manipulating the effective strength of 
the two-body interactions is based on a two-mode model for the condensate in the presence 
of an external trapping potential |10[ fl25| . This is discussed in appendix A. 



V. SUMMARY AND CONCLUSIONS 

We have shown that the atomic dressed states, which diagonalize the single body Hamil- 
tonian, are a useful basis for two-component condensates when the dressed state energy 
splitting is much larger than the mean- field energies. In this limit, the two-body interac- 



ts 



tions take on a simple form since those interactions, in the dressed basis, which change the 
total number of atoms in each of the two dressed states are highly suppressed. As such, 
the Hamiltonian in the dressed basis has the same form as the Hamiltonian in the original 
atomic basis with no external field coupling the two atomic states. However, in the dressed 
basis, the coupling constants for the two-body interactions are functions of the dressed state 
angle, a quantity which is experimentally controllable. Consequently, collective properties 
of the condensate which depend on the strength of the mean field interactions, such as the 
speed of sound, may be controlled by adiabatically varying the dressed state angle. The 
key requirement that hu c d ^> n U aa ,n Ubb,n U a b represents a purely technical challenge of 
achieving sufficiently high Rabi frequencies and detunings. 
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VII. APPENDIX A- ANISOTROPIC TRAP. 

Up to now we have limited ourselves to the case of no external trapping potential, 
Vi(r) = 0. However, another possibility for manipulating the strength of the two-body 
interactions involves using condensates in highly elongated cigar shaped traps. Let Z be the 
length of the condensate along the long axis, which we take to be in the z direction, and R 
the radius of the condensate. When the condition ^ » 1 and kR <C 1 is satisfied, then 
the propagation of sound is along the long axis of the trap and one may ignore changes in 
the radial direction. If one ignores the trapping potential in the z-direction and considers a 
radial harmonic potential with trapping frequencies u a and uj^ for the two-components, 

V t (p) = l -m^p 2 (41) 
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where p = ^/x 2 + y 2 is the radial distance from the trap axis, then one may express the field 
operators as 

Mr) = -±=J2a ap Mp)e ipZ/h ; (42a) 
V L p 

%(v) = ^-Y,<XhpMpY Pz/h \ (42b) 
V L p 

where if>i(p) are the ground state radial wave functions for the two states. Equations ( |42a| - 
[4 2 b| ) represent a two mode approximation with respect to the transverse coordinates. Phys- 
ically, this approximation consists of neglecting excitations above the ground state in the 
directions transverse to the axis of the condensate. 

When Tiuji 3> n Q Uki, the ipi(p) are the harmonic oscillator ground states 

(mu^l 2 _ mu)ip 2 j2n 



UP) = [-^) e- m ^' m . (43a) 

By substituting Eqs. ( |42a|442b| ) into Eqs. ( fTb~| - |6c| ) and using Eq. ( |43aj ), one obtains, after 
integration over the spatial coordinates, an effective one-dimensional Hamiltonian for the 
excitations in the axial direction which is given by, 

Hatom = E ( (J— + y + huJ aj a j ap a ap + - y + Hu b ] a\ p a bp + %& R (a\ p a ap + a\ p a bp ^j ] ; (44a) 



p 



H co ll — , \Uaad api C^ ap2 a a p 3 a, a p A + Ubb^b\)i a bp2 ab Ps a bp4 ^^ab a l l p 1 ( ^ b p2 0l ap3 ( ^bp4^ ■ (44b) 

P\+P2=P3+Pi 

One now has a rescaled Rabi frequency, Q R , and one dimensional mean-field interactions, 
Uij which are given by 



Cl R = Q R [ 2-npdpip a (p)ip b (p) = tt R 2v/UJaUJb ; (45a) 

U aa = U aa / 2irpdp\^ a (p)\ 4 = U aa -^; (45b) 

U bb = U bb f 2rcpdp\Mp)\ i = U bb ^; (45c) 
J 2-kTi 



U ab = U ab I 2rr pdp\Up)\ 2 \Mp)\ 2 = U ab ^-^^. (45d) 
J TrnuJa + u> b 

One sees that the Hamiltonian has the same form as Eq. ( pi| - Ph~D but for a one dimen- 
sional homogenous condensate. The harmonic oscillator energies appearing in H a tom can be 
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eliminated by a redefinition of the internal energy levels for \a) and |6). However, the one 
dimensional mean- field interactions, [7^, can now be adjusted by varying the trap frequen- 
cies. For example, if u)b/u) a = 2 and = U then Ubb/U aa = 2 and U a b/U aa = 4/3. This 
case is shown in Fig. 1(c). 

The system described by Eqs. ( |44a| - |44b| ) is different from the Thomas-Fermi limit con- 
sidered in [jrjl [ i~9fl . In the Thomas- Fermi limit, the radial wave function is a superposition 
of harmonic oscillator excited states, ip\ (p) = Y^ m c{p,U a a,Ubb,U a b)m,4'i,m( y p), where the 
i>i,m(p) is the m th excited state of the radial harmonic oscillator and the c(p, U aa , Ubb, U a b) m 
are functions of the chemical potential and coupling constants(see Eq. (|3])). Consequently, 
the resulting one dimensional form of H a tom will n ot be independent of the total number 
of atoms (which is determined by p) or the strength of the mean-field interactions in the 
Thomas- Fermi limit. 
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Fig 1. Plot of the speed of sound in units oiu — y n °^ l aa as a function of the dressed state 
angle 9. The dotted line is u c and the solid line is Ud- The dressed sate angle is restricted 
to the range < 6 < tt/2. (a) U bb /U aa = 3/2 and U ab /U aa = 2 (b) U bb /U aa = 1/2 and 
U ab /U aa = 2 (c) U bb /U aa = 2 and U ab /U aa = 4/3. 
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